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The results of the WA97 collaboration for strange particle production at mid-rapidity in Pb-Pb 
collisions at 158 GeV-A/c at CERN display a strong strangeness enhancement with system size 
at mid-rapidity which is dependent on the strangeness of the particle concerned, and saturates at 
values of participating nucleons greater than 120. These results are phenomenologically described by 
the mixed canonical ensemble, with canonical (exact) strangeness conservation involving all strange 
resonances, and grand canonical conservation of charge and baryon number. A detailed quantitative 
analysis shows that the data are well described by an equilibrium (7s = 1) hadron gas. 
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I. INTRODUCTION 



The WA97 collaboration § at CERN has shown that 
strange particle yields per wounded nucleon reach a sat- 
uration level for the most central Pb-Pb collisions (the 
Pb-Pb data is presented in 4 centrality bins) and show 
a pronounced increase when compared to a p-Be system. 
A mixed canonical description of these data has been 
published by Hamieh et al.[|[ ^ who showed that the 
density of strange particles derived from such a treat- 
ment depends on the size of the system consistent with 
the WA97 observations In addition Hamieh et al. ||] 
have shown that if reasonable values are used for the 
thermal parameters of each system the behavior of the 
yield (in the Pb-Pb systems) per wounded nucleon nor- 
malized to the p-Be yield is in agreement with the values 
obtained by the WA97 collaboration. The present work 
uses also he mixed canonical ensemble and presents 
a more thorough analysis in that the thermal parame- 
ters for each system are obtained by minimizing £ind 
a more explicit comparison between model and data is 
performed. The dependence of the parameters on the 
size of the system is presented, as well as the yields per 
wounded nucleon relative to p-Be. 



II. MIXED CANONICAL FORMALISM 

The mixed canonical partition function for a hadron 
jas conserving strangeness exactly is described by the 
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(1) 



following partition function j^: 

Z§=o = ^ £^ dcj^ exp p ' 

where S'„ = ^"^k^k^ ^ volume, and the sum 

is over all particles and resonances carrying strangeness 
n. For a particle of mass rrifc, with spin-isospin degen- 
eracy factor gk, carrying baryon number Bk and charge 
Qk with baryon chemical potential /is and charge chem- 
ical potential /ig, the one-particle partition function is 
expressed in the Boltzmann approximation as: 



Zt = — T^mti K2 



exp(Bfc/i_B + QfcAig) 



(2) 



As described in Ref. ji) the partition function for this 
ensemble may be rewritten in the following form which 
is well suited for numerical calculation: 



n -2n-3p 
2"l 



n— — 00 p— — 00 
Xln{x2)Ip{x3)I- 



2n—3p 



(a;i), 



where 



(3) 

(4) 
(5) 

and li are modified Bessel functions. 

The expression for the particle density, n^, may be 
obtained from the partition function Eq. (3) by follow- 
ing the standard method Q. For a particle i having 
strangeness s the result is 



Oi — \/ Si/ S-i 
Xi — 2-\/ SiS—i 



Z} 



\ ^ \ ^ p n -2n- 



3p — 3 



-00 p—~oo 



Xln{x2)Ipix3)I-2n-3p-sixi)- 



(6) 
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III. RESULTS 

In this section we explore whether canonical 
strangeness suppression at small volumes (compared to 
grand canonical equilibrium particle yields) is able to 
explain the enhancement of strange particle yields per 
wounded nucleon from small to large systems, as mea- 
sured by the WA97 collaboration at CERN||l|. The ther- 
mal model is well suited for to An integrated data and 
to the central rapidity region if a boost invariant plateau 
exists around at mid-rapidity We assume the valid- 
ity of the latter in order to analyze the particle yields. 
A thorough comparison of An particle numbers with the 
thermal model has been made in Ref. Q. As has been 
the case in other calculations ^ , our findings at mid- 
rapidity are found to be consistent with full strangeness 
equilibration (i.e. 7s = 1). This is in contrast to 47r 
integrated particle yields which are incompatible with 
full strangeness chemical equilibriun and deviate from 
it by several standard deviations. Experimentally, ratios 
of particle to anti-particle yields ||l5| at mid-rapidity ll6| 
have been compared to 47r integrated yields and shown 
to be in agreement for S+S collisions at 200 A GeV by 
the CERN NA35 collaboration @. In addition, _and 
perhaps of greater interest, the At: integrated ratio S/S 
measured by CERN NA49 for Pb-hPb collisions at 158 
A GeV has been shown to agree with the corresponding 
CERN WA97 mid-rapidity ratio fill. 

As the mixed canonical formalism has been derived 
here for the case of Boltzmann statistics, it is worth not- 
ing that the second term in a series expansion of the cor- 
rect quantum statistical distribution function for kaons 
gives a correction to the kaon numbers of the order of 
3% at a temperature of 150 MeV. The corrections to 
the Boltzmann distribution functions due to quantum 
mechanics are expected to affect the kaons more than 
any other strange particle, as they are the lightest of the 
strange particles. With errors of 3% and less, the use of 
Boltzmann statistics to describe the strange particles is 
justified. For all mixed canonical analyses the parameter 
fj,Q fit to zero, and has been subsequently removed as a 
free parameter. For a Att integrated Pb -I- Pb system, /ig 
is expected to be small and negative. 



A. Analysis A 



full quantum mechanical grand canonical particle number 
expression is evaluated for all non-strange particles. The 
parameters obtained in this way are shown in Table || as 
analysis A||l^. The for each system is good, and there 
are three (two) degrees of freedom for the p-f-Pb, and 
Pb+Pb systems (p-|-Be system). The large uncertainties 
in the radius parameter for the Pb-|-Pb systems shows 
that we have hit the grand canonical limit, and have no 
volume dependence, this is discussed briefly in a later 
section. 



B. Analysis B 

Analysis B has only three data points (^, § and =) 
for each fit. This analysis has been motivated by the rea- 
sons mentioned previously in the text. This analysis has 
no free parameters and shows reasonable agreement with 
analysis A. The large uncertainties of the fitted param- 
eters are unavoidable. Interestingly, fitting these ratios 
proves impossible for the p-|-Pb system. The hypoth- 
esis of 1^ of a special 'interaction volume' for strange 
particles in this system cannot be tested, as the parti- 
cle ratios would depend only on this interaction volume, 
as the regular volume cancels in the ratio of the particle 
multiplicities. 

C. Analysis C 

Analysis C determines the fits to the data using the 
grand canonical ensemble with quantum statistics. In 
this case, the volume dependence canels out in the par- 
ticle ratios and we have included /ig and /ig. This pro- 
cedure fails for the p-l-Be system. As may be expected 
for a system of this size, a canonical treatment [^] is 
required. In this analysis, the parameters fitted for Bin 4 
immediately catch the eye - with high uncertainty values 
even though < 1. This is especially puzzling when one 
notices that this is not the case in analyses A or B. The 
T and /is parameters agree with those obtained by Be- 
cattini et al. 0, for a grand canonical description of the 
47r-integrated data for central symmetric Pb^°* collisions 
at 158 A GeV by the CERN NA49 collaboration. The 
parameters in analyses A, B, and C agree within th each 
other one standard deviation. 



The main results of analysis A are based on fitting 
the ratios of the strange particle yields to the yield of 
negatives for each case in Table |. The negatives yield is 
used for each ratio due to the good statistics. In addition 
it provides a measure of the entropy of the collision. For 
negative pions, the quantum mechanical correction to the 
Boltzmann distribution function at 150 MeV is of the 
order of 30% . It is unacceptable to ignore this correction 
when attempting to describe the data. Considering that 
the non-strange particle yields predicted by the model are 
independent of the exact conservation of strangeness, the 



D. Thermal Model Parameters 

In order to fully reproduce the data obtained by the 
WA97 collaboration , it is necessary to determine a re- 
lationship between the parameters /is, the radius R and 
the temperature T obtained from the fits of analysis A 
and the number of wounded nucleons. The fits used 
are shown in Figs. ^ (temperature), and ^ (mb), as 
well as the parameters obtained from a full Boltzmann 
treatment J18[ . The volume paramaterisation is described 
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TABLE I: Parameters obtained for various analyses as described in the text. The temperature (T) and chemical potential (pb) 
are in MeV. The radii (R) are given in fm. If is greater than 50, the fit parameters are not shown. 





p+Be 


p+Pb 


Pb+Pb 








Binl 


Bin2 


Bin3 


Bin4 








Analysis A 








Ratios of Strange 


Particle to Negatives Multiplicities 




0.756 


1.85 


2.07 


0.427 


1.58 


1.45 


T 


162±4 


172±2 


161±4 


164±4 


166±4 


159±5 




111±7 


157±39 


201±25 


230±38 


228±28 


204±30 


R 


1.39±0.14 


1.16±0.39 


6.81±8.73 


9.88±10.9 


6.78±9.39 


10.4±6.8 




Analysis B 
Ratios of Particle Yields to that of their 


Anti-Particles 


x' 


~ 10"" 


60287 


~ 10~' 


~ 10~* 


~ 10"'' 


~ 10"^ 


T 


157±25 




171±17 


161±15 


153±35 


163±45 




106±31 




229±55 


222±46 


190±61 


220±30 


R 


1.50±0.83 




2.17±2.03 


6.06±8.92 


9.38±7.12 


8.08±6.83 








Analysis C 










Grand Canonical Fit 






x' 


78 


17.5 


1.53 


0.15 


0.70 


0.97 


T 




143±8 


173±20 


165±4 


166±5 


171±20 


fJ.B 




205±106 


282±96 


227±39 


203±45 


265±120 


fJ-Q 




-109±52 


-54±68 


0±50 


0±60 


-62±86 


MS 




105±63 


110±73 


54±25 


39±28 


94±97 



by equation |7|. The parameters obtained from the p+Pb 
data have not been considered because of their prediction 
of a volume smaller than that of the p+Be system. 



1. Temperature 

As may be seen in Fig. ^ it appears initially that the 
temperature increases with centrality. However, the last 
data point does not fit this hypothesis, and the varia- 
tion of chemical freeze-out temperature with number of 
participants has been assumed to be constant (T = 163 
MeV) . This is in agreement with the fit to the NA49 data 
by Becattini et al. 0| to central Pb-|-Pb data. This tem- 
perature value is slightly less than the one obtained by 
Hamieh et al. 1^. 

2. Radius 

In order to reproduce the WA97 data, knowledge of the 
variation of system size with average number of wounded 
nucleons is required. We have chosen a function of the 
form: 

< iVwound >=R^+b (7) 

with b < 0.5 adjusted to reproduce the data of the WA97 
collaboration (Figs. |^ and ||). This differs from the de- 
pendence of the radius on Apart assumed by Hamieh et 
al. @ where Apart ^ 1-3 - 1.7R^. 

The large uncertainty in the radius for the larger sys- 
tems is not surprising, as the ratios being fitted are ex- 
pected to show a small volume dependence, due only to 



175 I 1 , , 1 , 

Temp assuming full Boltzmann stats i — i — i 
Temp, assuming quantum stats, for non-strange particies ^--x--j 
Temp function usecf to reproijuce tire wA97 data 

170 - ] 



, i >t 

> 165 : ] 1 

2 ; ! ; 

2 1 i I 

3 160 I i 
CC 1 

i t 
Q. 

E T 



150 



145 I 1 ^ ^ 1 ^ ^ 1— 

50 100 150 200 250 300 350 

^'^wound'* 



FIG. 1: Variation of the chemical freeze-out temperature (T) 
of the system with number of wounded nucleons. 



canonical strangeness suppression. The effects of canon- 
ical strangeness suppression decrease with system size as 
one approaches the grand canonical limit where ratios of 
particle multiplicities have no volume dependence. 



3. Baryon Chemical Potential 



Figure |2| shows the variation of the baryon chemical po- 
tential with system size. This function increases rapidly 
before saturating. All the Pb-|-Pb bins are described 
by the saturation value of approximately 210 MeV. The 
obtained values of /zb in the central Pb-|-Pb bin are in 
agreement with that of Becattini et al. H . The /i used 
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FIG. 2: 
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Variation of the baryon chemical potential {hb) with 



for the Pb+Pb bins by Hamieh et al. is just outside 
one standard deviation of the values obtained during this 
analysis. For the p+Be system the value of /is used in 
Ref. (150 MeV) is much larger than that obtained in 
this analysis (111 ± 7 MeV) . 

The exact variation of /xb, T and R in the range be- 
tween the p-Be and Pb-Pb may differ from what is shown 
here and new data from CERN NA57 in this range are 
eagerly awaited. The figures mentioned above (111 and 
^ also show the variation of the parameters T, and /is 
with centrality in the case where all particle multiplic- 
ities are assumed to be Boltzmann. These parameters 
are seen to be in agreement with the parameters from 
Analysis A. The value of the parameter obtained for 
the purely Boltzmann fits were of order 1 for all systems 
considered. 



E. Predicted Evolution of Particle Yields 

Figs. H and |4| show the ability of the model to repro- 
duce the data using the functions in Figs ^ and and 
equation |^ to predict the variation in chemical freeze-out 
temperature T, baryon chemical potential fiB, and ra- 
dius R with average number of wounded nucleons. As is 
clear, the agreement of the model with experimental data 
is good. The exact shape of the thermal model predic- 
tions shown in Figure ^ and ^ show a large dependence 
on the exact relationship between number of wounded 
nucleons and system size, and a weaker dependence on 
the functional fit to the other two parameters. This is 
rather unfortunate as the thermal model R values ob- 
tained have uncertainties equal to their magnitude, so 



this relationship is not constrained by the model. 
IV. CONCLUSIONS 

A thermal model conserving baryon number and 
charge on average, and strangeness exactly has been used 
to describe particle yields from heavy ion collisions. The 
formalism includes all strange particles. The model has 
been appHed to the CERN WA97 data, and is shown to 
be able to reproduce the data for all centrality classes 
(Figs. |and|). The accuracy of the model parameters 
(temperature T, radius i?, and baryon chemical poten- 
tial ^s) is restricted by the lack of 47r integrated particle 
yields for multi-strange particles. In order to reproduce 
WA97 graphical presentation of their data accurately, 
some knowledge of the variation of T, R, and /is was 
required. A functional form of the variation of these pa- 
rameters with the number of wounded nucleons has been 
determined and is presented in Figs and ^, and equation 

0- 

The parameters obtained in this model for the cen- 
tral Pb-|-Pb bins are in agreement with the An thermal 
model application of Becattini et al. 0. This lends some 
weight to the accuracy of the model when fitting particle 
ratios in a limited kinematic region. This grand canoni- 
cal model includes the factor 75 ^ 1 to predict the yields 
of strange particles. To differentiate between canonical 
strangeness suppression, and suppression of strange par- 
ticles by anomalous phase space occupancy fl^ , the yield 
of the (f) particle could be used. Yields of this particle are 
not sensitive to canonical strangeness suppression, but as 
it contains an s and an s quark, these yields will be sen- 
sitive to 75 ^ 1. 

The enhancement of strange particles measured by 
CERN WA97 has been considered a signal for QGP for- 
mation. At first glance, the ability of a full equilibrium 
thermal model to reproduce the data suggests the ex- 
istence of a deconfined state, since equilibrium strange 
particle yields have been proposed as a possible signal 
for deconfinement ||l^. A deconfined phase is not, how- 
ever, expected to be formed in the p-fBe system, as it 
is only expected in large dense systems. The ability of 
the model to reproduce the p-f Be data suggests that it is 
possible for strange particles to reach equilibrium yields 
by hadronic interactions alone. 
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